Robustness of the nodal rf-wave spectrum to strongly fluctuating competing order 
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We resolve an existing controversy between, on the one hand, convincing evidence for the existence 
of competing order in underdoped cuprates, and, on the other hand, spectroscopic data consistent 
with a seemingly homogeneous d-wave superconductor in the very same compounds. Specifically, 
we show how short-range fluctuations of the competing order essentially restore the nodal d-wave 
spectrum from the qualitatively distinct folded dispersion resulting from homogeneous coexisting 
phases. The signatures of the fluctuating competing order can be found mainly in a splitting of 
the antinodal quasi-particles and, depending of the strength of the competing order, also in small 
induced nodal gaps as found in recent experiments on underdoped La2- x Sr x Cu04. 
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Many recent experiments point to the prominence 
of competing phases in underdoped cuprate supercon- 
ductors, and challenge existing theories for the pseu- 
dogap phase, [l- 12 1 Evidence has been found for stripe 



14| and YBa 2 Cu 3 6 



phases in the La2Cu04[l 

families, [III, ll5Ml7| checkerboard and nematic phases 
in the Bi-based cuprates. (l8T - l2Ci | loop-current order in 
YBa2Cu306+a;,[2l| and spin glass phases in most of 
the highly underdoped cuprates. }22l424| All of these 



phases are expected to have clear spectroscopic signa- 
tures which, in many cases, involve a Fermi surface re- 
construction; puzzlingly, however, angle resolved pho- 
toemission spectroscopy (ARPES), which has become 
a powerful probe of the electronic structure of the 
high temperature superconductors . finds no clear evi- 
dence of any reconstruction. [2514281 ] even for materials 
like La2- x Sr x Cu04 (LSCO) near x = 1/8 where signa- 
tures of reconstruction due to stripes are expected to be 
maximal. [2^, [30[ A few recent studies report antinodal 
(AN) particle-hole symmetry breaking in the pseudogap 
phase of Pbo.55Bii, 55 Sri. 6 Lao.4Cu0 6+ 5 (Bi2201),[iB, 111 
thus confirming the presence of competing order, but the 
underlying low-temperature dispersion of the electronic 
quasi-particles in the nodal region remains remarkably 
similar to that of a well-known c?-wave superconductor. 
This is still true even in the insulating spin-glass regime 
at doping levels so low that superconductivity is not yet 
present. [27j In many cases, the competing order appears 
to be fluctuating rather than static, 0, |22[ and in this 
paper we demonstrate that many aspects of the ARPES 
spectrum, such as the robustness of the nodal d-wave 
spectrum to strong competing order, can be understood 
if this fact is accounted for. 

Motivated by the widespread observation of short- 
range slow antiferromagnetic (AF) fluctuations in under- 
doped cuprates, we focus on the competition between 
AF and d-wave superconducting (dSC) order, but we ex- 
pect our main results may apply also to other candidates 
for the competing order. Specifically, we have performed 
Monte Carlo (MC) simulations of a two-dimensional dSC 
in which thermal fluctuations of both the AF and dSC 



order parameters are included. Our main finding is that 
highly inhomogcncous phases, which are widely reported 
experimentally, qualitatively alter the low-energy quasi- 
particlc dispersion from the homogeneous or weakly dis- 
ordered situation. We find that: (1) unless the compet- 
ing order is very strong, the spectrum is indistinguishable 
from that of the pure dSC, (2) even for strong competing 
order the nodal electronic spectrum is very similar to that 
of the pure superconductor, and shows no sign of Fermi 
surface reconstruction expected from AF order, and (3) 
when the competing order is strong, the AN gap has a 
distinct energy scale from the nodal region. These re- 
sults bridge the seemingly contradictory findings of sim- 
ple momentum-space behavior reported by ARPES ver- 
sus complex inhomogeneous real-space behavior reported 
by local probes. 

Our main results are illustrated by Fig. [TJ where spec- 
tra are shown for a pure dSC and for a dSC with strong 
competing magnetism (dSC+AF). The top panels (a-d) 
show snapshots, taken from the MC simulation, of the 
dSC amplitude A-d{ r i) (b,d) and of the local AF moment 
mg(rj) (a,c). The bottom panels (f,g) show the spectral 
function A(k, lu) at k-points taken along the nodal and 
AN lines shown in [TJe) . The spectral function shown in 
Fig.QIf) [Fig.QJg)] arises from sampling of configurations 
similar to (a,b) [(c,d)]. For comparison, in Fig. [ljh) we 
also show spectra for the simplified case when dSC and 
AF order coexist homogeneously In this case, we use 
the root-mean-square values of the AF and dSC order 
parameters fromHJg). 

The similarity between the nodal spectra for the pure 
dSC [Fig. Hp)] and dSC+AF [Fig. QJg)] is remarkable, 
especially in comparison to the uniformly coexisting case 
[Fig. Hp 1 )]. I 11 the uniform case, the nodal cut exhibits 
the well-known band back-folding that results in the for- 
mation of closed Fermi surface pockets around the AF 
Brillouin zone boundary. The back-folding is due to the 
large AF-induced band gap that appears above the Fermi 
energy. It is striking that the band gap is completely 
absent inQJg); instead, the dispersion is almost indistin- 
guishable from that of the pure dSC despite the presence 



2 




(5) °° l(d) ' " 




0.3 







It 


-0.3 




0.6 




0.3 


-31 








(e) 



J AN 



Uniform SC+AF 




-2-1 1 -2-1 1 2-2-1 1 -2-1 1 2 

CO CO CO 0) 



-2-1 1 -2-1 1 2 

CO CO 



FIG. 1: (Color online) (a-d) Snapshots of MC simulations 
for (a,b) a pure dSC and (c,d) a dSC with competing AF 
order at reduced temperature r = T/T c — 0.3. (a,c) The lo- 
cal AF moment mQ(ri) and (b,d) dSC order parameter are 
shown for the final step in the MC simulation, (f-h) The av- 
eraged spectral function is shown along nodal (N) and antin- 
odal (AN) cuts illustrated in (e). The spectra are for (f) dSC 
only, (g) dSC+AF, and (h) a homogeneously coexisting dSC 
and AF. Linewidths are due to thermal fluctuations, except 
in (h) where they are artificial. In all cases, the filling is 
n ~ 0.88. Parameters are V = —1.0 and (a,b) [7 = 0, for 
which T c = 0.20 and (c,d) U = 3.5, for which T c = 0.06. The 
Fermi energy is ef — throughout. In (h) we set Ad = 0.35 
and UrriQ — 3.5 x 0.2 based on average values in (g). 



of sizable magnetic moments on nearly all sites evident 
from Fig.QJc). The only signature of the competing mag- 
netism is the narrow nodal gap at the Fermi energy in 
[ljg) similar to what has been recently found in ARPES 
studies of underdopcd LSCO.[32| We stress that in our 
model, the nodal gap appears only for large values of 
U, close to the superconductor- insulator transition; at 
smaller values of U, the nodal dispersion remains indis- 
tinguishable from the pure dSC case. 

In contrast to the nodal spectrum, the AN spectrum 
in (g) resembles that of the uniform dSC+AF case in 
(h): the single quasiparticle peak at oj < ep in (f) is 
split into a pair of peaks in (g) and (h), one of which is 
upward-dispersing and the other of which is downward- 
dispersing. While the splitting is obvious in (h), the 
peaks are not easily resolved in (g) because the line 
shapes are thermally broadened, and the lower intensity 
peak appears as a shoulder. 



The results shown in Fig. QJg) are qualitatively similar 
to what has been reported at low T for ARPES exper- 
iments on Bi2201, lfl 3l| namely, there is a band split- 
ting at the antinode but no sign of Fermi surface recon- 
struction near the node. We emphasize that, while it 
is possible to model the AN spectrum with a uniform 
"finite-t7" spin or charge density wave, as in Ref. |3lj, it 
is much more difficult to model the spectrum along the 
entire Fermi surface that way. Indeed, Norman et al. j33| 
have argued against finite-^ models of the pseudogap for 
this reason. Figure QJg) shows that fluctuations qualita- 
tively alter the spectrum in a way that is consistent with 
experiments. 

Figurc[T]is based on MC simulations on an L x L lattice 
of electrons coupled to two thermally fluctuating classical 
fields. The real field hi couples to the local magnetiza- 
tion at site i while the complex field dij couples to Cooper 
pairs along bonds connecting nearest-neighbor sites i and 
j. Purely superconducting models (with hi = 0) have 
been widely used to study phase fluctuations, 34hl38j but 
relatively little work has focussed on the more compli- 
cated problem of competing order. The partition 
function is 



Z = J D[hdd*]exp[-/3n(h,d)}, 



(1) 



where Sl(h,d) = -Tin Tr exp(-pH) 

over electronic degrees of freedom, J D 
dimensional integral over the fields, and 



, Tr ... is a trace 
hdd*} ... is a 5L 2 - 
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In Eq. Uj are the hopping matrix elements between 
nearest (Uj = —t= — 1 ; i is the unit of energy throughout 
this work) and next-nearest (tij = t' = 0.4) neighbors. 
To reduce the number of integration variables, we impose 
the singlet constraint dij = dji. U and V control the size 
of the AF and SC fluctuations respectively; at T = 0, 
the saddle point approximation is exact and gives hi = 
Um(Yi) where m(r<) = and d tj = V{cj±Ci\). 

Throughout this work, we fix the pairing interaction to be 
V = —1.0 and vary the Coulomb repulsion U. The dSC 
order parameter is Ad(r ? ) = ^2j{— ^) Vi ~ Vl {cj^c^), where 
j is summed over nearest neighbors of i, and the AF mo- 
ment is m(j(r,;) = (— l) Xi+yi m(Yi). The dSC transition 
occurs at the temperature T c where the pair correlation 
A d (r t )A* d {r t + R) = 0, with R = \{L, L). 

We use a Metropolis algorithm to evaluate the integral 
in Eq. ([T]). To reduce noise, all spectra are averaged over 
10-20 separate MC simulations; a typical simulation con- 
sists of 2 x 10 4 sweeps, where one sweep has L 2 steps, 
in which each field is updated. The MC calculations are 
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FIG. 2: (Color online) Results of MC simulations at T = 0.02 
for three different U values. Results are for (a-d) U = 3.3, 
(e-h) U = 3.5, and (i-1) U = 3.8. First and second columns 
contain MC snapshots of mg(ri) and A<j(ri) respectively; the 
third column contains the LDOS corresponding to the config- 
uration in the first two columns. The final column shows the 
MC-averaged spectral function .A(k, ep). 



computationally demanding because H must be diago- 
nalized at each step; in similar calculations [39f systems 
were limited to L < 12 sites (see, however, [37| for an 
alternative approach). Here, we use first order perturba- 
tion theory to estimate updates to Sl(h, d) at each MC 
step, with a full diagonalization of H at the beginning 
of each MC sweep. We have checked that this does not 
change the MC-generated distribution of hi and dij val- 
ues, and can study systems up to L = 32 by this method. 

Typical mg and Ad configurations sampled by the 
MC simulations shown in Fig. [2] illustrate how the model 
evolves with increasing U. For U = 3.3 [Fig. GJa)], there 
are large regions where TOg(rj) is small, and small pock- 
ets of short-lived AF order. The dSC order parameter 
also fluctuates [Fig. |2Jb)], and there is a clear spatial 
anticorrelation between A^r^) and mg(r i ). In^c), we 
have plotted the local density of states (LDOS) at the 
Fermi energy, calculated for the particular configuration 
of m Q {vi) and A d (r l ) shown in^a) and^b). The LDOS 
is reduced inside the AF pockets, and is largest in re- 
gions where the AF moments are small. Finally, in[2][d), 
we show the MC-averaged spectral function at ef- It is 
apparent that the nodal quasiparticles are largely unaf- 
fected by the fluctuations, so that the low energy spectral 
weight is concentrated at the four nodal points. 

The AF fluctuations are larger when U = 3.5 
[Fig. Efe)], and A^ is suppressed in a significant fraction 
of the sample, although the maximum value of A^ in[2jf) 
is about the same as in[2jb). The tendency for the low 
energy LDOS, shown in^g), to be concentrated along 
AF domain walls is more pronounced than in[2Jc), and 
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FIG. 3: (Color online) Spectral function A(k, u) along the 
Fermi surface. Curves (offset for clarity) are taken at k-points 
between the nodal point (top curve) to the AN point (bottom 
curve). Panels (a)-(e) correspond to different U or reduced 
temperature r = T/T c , and are arranged in order of increasing 
average magnetic moment. T c values are shown in (f); root- 
mean-square averaged (g) m and (h) Ad is shown at r = 0.3 
(solid symbols) and r = 0.5 for U = 3.5 (open symbols). 



becomes even more pronounced when U = 3.8, Fig.[2][k), 
where there is now phase separation. In this case, mo- 
bile holes lie almost entirely at the boundaries between 
AF domains; there is some residual pairing of the holes 
[Fig. HUj)], but there is no long-range phase coherence 
and the system is non-superconducting. 

The spectral function at ef for U = 3.8 shown in 
Fig. (2^1) reveals very little of the complex real-space 
structure that emerges as U increases. Indeed, the main 
change is that the nodal points in Fig. [D[d) evolve into 
"Fermi arcs" with increasing U. That AF fluctuations 
contribute to arc formation was suggested previously in 
Ref. 0; here we add the observation that this persists 
even into the nonsupcrconducting state. Indeed, it is re- 
markable that one recovers the underlying "bare" Fermi 
surface purely from states along the spaghetti-like do- 
main walls. This situation is reminiscent of the low- 
energy spectral features studied in the case of disor- 
dered stripes, 42- 44j which, however, assumed perfect or- 
der along one spatial direction leading to a characteristic 
Fermi surface reconstruction. 

Next, we show in Fig. [3] the progressive evolution of 
A(k, u>) as the magnetic fluctuations are increased. The 
spectra are taken at momenta along the Fermi surface be- 
tween the nodal and the AN points. To obtain good mo- 
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FIG. 4: (Color online) Spectral gap along the Fermi surface 
for different U and reduced temperature r. The angle 6 is 
measured along the Fermi surface, with 6 = corresponding 
to the antinode and 9 = 45° to the node. Inset: Decompo- 
sition of the antinodal spectrum into gaussians for U = 3.5, 
t = 0.5. 



mentum resolution, wc calculated the spectrum at points 
that are interpolated between the L 2 k-points of the orig- 
inal simulation as follows: after the MC simulations were 
completed, we recalculated the spectrum with a com- 
plex boundary condition using the previously-generated 
sequence of hi and dij values. The complex phases were 
tuned to shift the k-values to the Fermi surface (a differ- 
ent phase is required for each k-point). In panels (a-c) 
A(k, w) is shown for increasing U at the same reduced 
temperature r = T/T c = 0.3, and the root-mean-squarc 
values of and m at this r are plotted in Etg) and 
Et n ): respectively. InEtd), r = 0.5, which increases m 
and lowers Ad, as shown by the open symbols inEtg) and 
Eth) . Finally, in [3je) , spectra are shown for the phase- 
separated system at low T. 

In agreement with the main conclusions from Fig. [TJ 
Fig. EJb) is essentially indistinguishable from Fig. Eta), 
even though T c is reduced by a factor of two by AF fluc- 
tuations [Fig. 024)]. Small differences from the <i-wave 
spectrum only emerge in Et c ) , by which point T c is one 
third of its value in[3ta). In Etc), a small gap appears 
at the node and the AN peak below ep splits into two. 
These features are more pronounced in Etd) , where the 
AF fluctuations are larger. Finally, in Etc), there is a 
reorganization of the electronic structure into lower and 
upper magnetic bands. It is the small residual spectral 
weight near ep, coming from the AF domain walls, which 
generates the Fermi surface in Fig. [2]f 1) . 

As mentioned above, a recent comprehensive ARPES 
study of LSCO has discovered that the excitation spec- 
trum is gapped along the entire underlying Fermi sur- 
face, and hence displays a nodal gap, in the highly un- 
dcrdopcd regime. [32J At doping levels where the nodal 



gap is observed, a spin glass is known to coexist with 
superconductivity. [22f and in Fig. [4] we show that AF 
fluctuations can indeed generate a nodal gap. 

Figure[4]shows how the gap along the Fermi surface de- 
pends on both the magnitude and the correlation length 
of the AF fluctuations. To be consistent with experi- 
ments, we define the gap at a point k on the Fermi surface 
as the position of the highest-lying peak in A(k, to) for 
ui < ef- When U < 3.3, the gap along the Fermi surface 
has a dSC-like structure, and is characterized by a single 
energy scale. When AF correlations are stronger (see e.g. 
the U = 3.5, r = 0.3 curve) the AN gap remains set by 
the dSC scale, and a second energy scale emerges in the 
form of a small gap at the node, similar to what was found 
by Razzoli et al. [321 ] In Fig. [4] the nodal gap grows with 
the AF magnitude, but vanishes at high temperature 
when the AF correlation length is short and line broaden- 
ing wipes it out. A third energy scale emerges at the AN 
when AF correlations are strong (U = 3.5 and t > 0.5). 
In this case a large gap, with an energy scale distinct 
from that near the node, develops. Similar crossovers 
between a superconducting nodal gap and nonsupercon- 
ducting AN gap have been widely seen in ARPES exper- 
iments on many underdoped cuprates;[l[ 0, [|1 0, E3; 12 1 
that this is not reported in LSCO suggests that AF cor- 
relations are not strong enough for a clear AN feature to 
develop. 

We have found that, even when they cannot be resolved 
by eye, two energy scales are always present in the AN 
spectrum when U is large (Fig. [4] inset). The inner peak 
has a d-w&ve k-space structure, while the outer peak is 
tied to the AF fluctuations, although it is not simply 
related to the magnetic energy scale Um ~ 0.7. To con- 
firm AF correlations as the source of the nodal gap in 
LSCO, our calculations suggest that one should look for 
a second energy scale, particularly in the AN spectrum. 
This would distinguish spin glass physics from models 
(such as d + id pairing [32 1) in which the gaps are added 
in quadrature. 

In conclusion, we have shown that the overall spectrum 
measured by ARPES experiments is consistent with sim- 
ple models of competing order if one accounts for the 
highly inhomogeneous nature of the competing phases. 
Such models explain the origin of the robust electronic 
states in momentum space arising from complex real- 
space environments with substantial amplitudes of the 
fluctuating competing order. 
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